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CHEMICAL KINETIC MODELING OF PROPANE OXIDATION BEHIND SHOCK WAVES 


Allen G. McLain and Casimir J. Jachimowski 
Langley Research Center 


SUMMARY 

The stoichiometric combustion of propane behind incident shock waves was 
studied experimentally and analytically over a temperature range from 1700 K to 
2600 K and a pressure range from 1.2 to 1.9 atm. Measurements of the concentra- 
tions of carbon monoxide [CCf] and carbon dioxide [COg] and the product of the 
oxygen atom and carbon dioxide concentrations QO] [CO] were made after passage of 
the incident shock wave. A kinetic mechanism was developed which when used in a 
computer program for a flowing, reacting gas behind an incident shock wave pre- 
dicted experimentally measured results ( QCO] , [CO2] , CQ] dCCQ, and reaction times) 
quite well. Ignition delay times from the literature were also predicted quite 
well. The kinetic mechanism consisted of 59 individual kinetic steps. 


INTRODUCTION 

The ability to model successfully the kinetics of hydrocarbon oxidation, 
particularly for high-molecular-weight hydrocarbons, is important in the develop- 
ment of reliable mathematical models of combustion systems, such as aircraft gas 
turbine combustors. Even though much work has been done in the area of combus- 
tion chemistry, very few experimentally validated chemical kinetic models are 
available. However, as a result of such work, guidelines for the assembly of 
hydrocarbon combustion models do exist. 

Recent experimental studies of the oxidation of high-molecular-weight hydro- 
carbons have indicated that these long-chain molecules decompose rapidly during 
early stages of combustion into hydrocarbon molecules and fragments with lower 
molecular weights. For instance, reference 1 indicated that larger hydrocarbons 
(iso-octane and n-heptane) degraded rapidly into ethylene, carbonyls, C3 ole- 
fins, and acetylene. Interestingly, this rapid decomposition was also noted in 
a study (ref. 2) of the oxidation of ethane, the next hydrocarbon above methane 
in the alkane series. After scrutiny of hydrocarbon combustion literature, 

Cooke and Williams suggested in reference 3 that rapid decomposition of the par- 
ent fuel followed by oxidation of the low-molecular-weight hydrocarbon molecules 
and fragments best described mechanistically the combustion behavior of long- 
chain alkane hydrocarbons. Because of these observations, a kinetic rate mecha- 
nism which can explain the combustion behavior of a high-molecular-weight alkane 
(propane, for instance) must incorporate kinetic steps observed for simpler 
hydrocarbons (CH4, C2H2, C2H4, and CgHg). In reference 4, methane oxidation 
behind shock waves was studied and a kinetic rate mechanism was developed which, 
when used in an analytical model for a flowing, reacting gas behind a shock 
wave, gave results that agreed well with experimental data. In a subsequent 
study (ref. 5), both acetylene and ethylene oxidation behind shock waves was 
studied over the temperature range from 1700 to 2600 K, and a kinetic rate 



mechanism was developed. This rate mechanism, incorporating kinetic steps from 
methane combustion, gave analytical results which again compared well with experi- 
mental data. The logical extension of this work is the study of the combustion 
of propane and the development and verification of a kinetic rate mechanism in 
which the propane decomposition reactions are coupled with the reactions that 
described oxidation of fragments with lower molecular weights, such as the 
methyl radical, acetylene, and ethylene. 

The decomposition and ignition characteristics of propane have been studied 
extensively in references 6 to 10. In reference 11, Chinitz and Baurer pre- 
sented a lengthy kinetic rate mechanism for propane in which many of the rate 
constants were calculated from the kinetic theory of gases and empirical rela- 
tions for estimating the Arrhenius preexponential factor and the activation 
energies. A number of the elementary reactions in this mechanism appear in the 
methane, acetylene, and ethylene kinetic mechanisms mentioned previously. How- 
ever, this lengthy mechanism assumed the formation of a number of oxygenated 
intermediates (not observed experimentally) and also relied on radical reactions 
as the primary path for destruction of the propane. 

The purpose of this paper is to present the results of the stoichiometric 
oxidation of propane behind shock waves at temperatures between 1700 K to 
2600 K and pressures between 1.2 to 1.9 atm. The progress of the oxidation 
reactions behind the shock waves was monitored by measuring the emitted radia- 
tion from carbon monoxide, carbon dioxide, and the chemiluminescent reaction, 

0 + CO -*■ CO 2 + hv. The measured intensity of the radiation was used to calcu- 
late concentrations of carbon monoxide and carbon dioxide and the product of the 
oxygen atom and carbon monoxide concentrations. A kinetic rate mechanism was 
assembled, and with the use of the analytical model of reference 12 for a flow- 
ing, reacting gas behind incident shock waves, a parametric study of the kinetic 
mechanism was made. Predicted and measured concentrations and reaction times 
were compared during the study resulting in the evolution of a kinetic mechanism 
for propane. Ignition delay times were predicted by the developed kinetic mech- 
anism and were compared with experimental ignition delay times from reference 6. 


SYMBOLS 

[CO], [CO 2 J molar concentrations of CO and CO 2 , mol/cm3 

[0] [CO] concentration product of atomic oxygen and carbon monoxide, (mol/cm3)2 

voltage measured with specific detectors and filters of wavelength A, 
V 

hv radiation emitted during chemiluminescent reaction, 0 + CO CO 2 + hv 

k rate coefficient, sec - ^ for unimolecular reactions, cm3/mol-sec for 

bimolecular reactions, cm^/mol^-sec for termolecular reactions 

P 
T 


2 


pressure, atm (1 atm = 101. 3 kPa) 
absolute temperature, K 



T m temperature predicted by kinetic model, K 

t-^ ignition delay time, ysec 

t m reaction time, ysec 

Abbreviations: 

IR infrared 

UV ultraviolet 


EXPERIMENTAL APPARATUS AND MEASUREMENTS 

All experiments were carried out behind incident shock waves in a stainless- 
steel shock tube with an inside diameter of 8.9 cm. Approximately half of the 
tests were carried out in a 580-cm-long test section with an observation station 
located 565 cm from the diaphragm location (station (T) in fig. 1), and the 
other half were carried out in a 671-cm-long test section with an observation 
station located 625 cm from the diaphragm location (station (2) in fig. 1). 

This test section and observation stations were used in the studies of refer- 
ences 4 and 5. 

Figures 2(a) and (b) are graphical illustrations of the 565-cm and 625-cm 
observation stations, respectively. Infrared monitoring equipment at the 625-cm 
observation station (fig. 2(b)) consisted of two indium antimonide infrared 
detectors mounted on opposite sides of the shock tube. At 90° from this loca- 
tion, a photomultiplier and a pressure or heat-transfer gage were mounted on 
opposite sides of the shock tube. The infrared detectors and the photomulti- 
plier were rigidly mounted to the shock tube with an adapter which contained a 
calcium fluoride window, collimating slits, and an interference filter of the 
appropriate wavelength for observing emission from CO, CO 2 , or the chemilumines- 
cent reaction, 0 + CO CO 2 + hv. The interference filter used to monitor CO 2 
radiation was centered at 4.3 yin with a bandwidth at half-peak of 0.15 ym, while 
an interference filter centered at 5.0 ym with a bandwidth at half-peak of 
0.15 ym was used to monitor radiation from CO. The interference filter used 
with the photomultiplier (having an S-20 response) for monitoring the radiation 
from the reaction, 0 + CO ** CO 2 + hv, was centered at 3700 A with a bandwidth 
at half-peak of 65 A. The fourth location at this observation station was 
equipped with either a dynamic pressure gage or a heat-transfer gage which could 
be used to trigger counters or to indicate shock arrival at the observation 
station. 

In the 580-cm-long test section (fig. 2(a)), only two measurements (emis- 
sion from the reaction, 0 + CO CO 2 + hv, and emission from either CO or CO 2 ) 
could be made behind the incident shock wave for each run. Thus, two runs were 
required to collect [ 6 ] [CO^ , QCO} , and [CO 2 ]] data. In addition, a monochromator 

was used to collect the radiation at 3700 A for the chemiluminescent reaction, 

0 + CO CO 2 + hv. The monochromator was equipped with a grating blazed at 
3000 A and a photomultiplier with an S-20 response. The slits were set to allow 



a half-peak bandwidth of 32 A at the 3700-A setting. The windows were calcium 
fluoride, and a quartz lens with a focal length of 80 mm was used to pick up the 
radiation from the center of the tube and focus it on the entrance slit of the 
monochromator. Slits placed between the calcium fluoride windows and the mono- 
chromator reduced the field of view to provide less than 2 usee of time resolu- 
tion in flowing gas. The infrared equipment used to monitor the emission from 
CO and CO2 was the same as for the 625-cm observation station. 

The following calibration procedure was used at both the 565-cm and 625-cm 
observation stations. The infrared detection system was calibrated by shock 
heating a mixture of CO or CO2 diluted with argon and recording the intensity of 
the emitted radiation as a voltage output on an oscilloscope. A quantitative 
relationship was established between the intensity of the -radiation measured 
from the CO or CO2 and the concentration of CO or CO2 in the test section over 
a range of temperature from 1700 K to 2600 K. 

The voltage output and the calibration curves for the C0 2 and CO infrared 
detectors are presented in figures 3 and 4, respectively. Figures 3(a) and 4(a) 
illustrate the variation in voltage output from the infrared detectors equipped 
with the 4.3-Mm and 5.0-ym filters for different mole fractions of CO2 or CO in 
the test gas mixtures. Figures 3(b) and 4(b) illustrate the empirical quantita- 
tive relationships developed between [CO2] and [C0j and the measured voltages 
and temperatures. Since the CO2 emission band overlaps the CO emission at 
5.0 pm, a correction was necessary to obtain the true CO concentrations. To 
make this correction, a correlation between the temperature and the ratio of the 
voltages measured at 5.0 pm and 4.3 Pm was determined for various C02~Ar mix- 
tures. This correlation is shown in figure 5. 

For the reaction, 0 + CO -*■ CO2 + hv, the detection system was calibrated 
by shock heating a mixture of H2-02-C0-C02~Ar over the same temperature range 
and recording the chemiluminescent emission. Details of this calibration pro- 
cedure are given in references 4 and 5- 

The output voltage with respect to time was displayed on an oscilloscope 
and photographed. The time constant for the infrared recording system was about 
3 psec. For the ultraviolet recording system, the time constant was less than 
2 psec. The speed of the shock wave was measured with a raster system upon 
which timing marks were superimposed at 10- or 50-psec intervals. On this same 
raster, output voltage deflections from heat-transfer gages located 30.48 cm 
apart indicated shock passage. In addition, the heat-transfer gage output at 
station (T) was used to start a counter operating in tenths of microseconds. 

The heat-transfer gage or dynamic pressure transducer at observation station (2) 
yielded a voltage output which was used to stop the counter. The time for the 
shock wave to pass from station (T) to the observation station (2) was used to 
calculate the shock velocity at the instant of arrival of the shock wave at the 
observation location. The times from the raster and the counter compared well. 

The test gas mixture was prepared by the method of partial pressures from 
commercially prepared propane (99-5 percent), oxygen (99.7 percent), and argon 
(99.995 percent) without further purification. Samples of the mixture were ana- 
lyzed with a commercially available hydrocarbon analyzer for the percentage of 
propane in the mixture which was found to vary not more than +0 . 1 percent from 


4 
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the desired value. The initial test gas pressure was 40 torr (1 torr = 133-3 Pa) 
for all tests. 


EXPERIMENTAL RESULTS AND DATA ANALYSIS 

An example of the experimental data obtained by shock heating a mixture 
(1 -percent CgHg, 5-percent O2, and 94-percent Ar by volume) is shown in fig- 
ure 6. The traces represent the oscilloscopic display of voltage output of the 
\ [CO] [0] , [CO] , and [002^ optical detection equipment as changes occurred with 

time. These traces are typical of the mixture being studied. All the [CO] [0] 
y chemiluminescent ultraviolet emission profiles at 3700 A exhibited a pronounced 

spike a short time after shock arrival followed by a slow decrease. Similar fea- 
tures were observed in references 4 and 5 during the combustion of methane, 
acetylene, and ethylene behind shock waves. The infrared emission profiles for 
CO2 at 4.3 pm and CO at 5.0 pm always exhibited a rapid increase after shock 
arrival followed by a knee and then a region of very slow increase. The knee of 
each infrared emission profile occurred slightly after the 3700-A emission 
spike. 

The 3700-A emission profile was analyzed by the method outlined in refer- 
ence 4. The entire postspike emission was assumed to be from the chemilumines- 
cent reaction, 0 + CO -*■ CO2 + hv. The maximum in this reaction's emission was 
assumed to occur at the same time as the maximum for the emission spike. By 
extrapolation of the postspike emission into the spike region, a maximum for the 
0 + CO -*■ CO2 + hv reaction profile was obtained. The voltage at this maximum 
was used to calculate the concentration product [0] [CO] . The [0] [CO] values cal- 
culated by this method and the reaction times t m defined as the interval (lab- 
oratory time) between shock arrival and the emission maximum are presented in 
table I . The uncertainty in calculating the [0] [CO] product in this manner was 
estimated to be 5 to 10 percent from the analysis in reference 4. The tempera- 
tures T in table I correspond to the temperatures immediately behind the shock 
before any reactions in the system could occur. These temperatures were calcu- 
lated from the computer program of reference 12 given the experimentally 
obtained shock velocity and assuming frozen shock conditions. Although the tem- 
perature at the time of the maximum in the [0] [CO] profile is greater than at 
the shock front, a temperature correction was not deemed necessary since the dif- 
ference in temperature would reduce the [o] [CO] product by less than 3 percent 
which is within the experimental scatter of the data. 

For reduction of the infrared data, however, the temperature at the maximum 
of the [0] [CO] profile must be used. The [CO] and [CO2] concentrations varied 
by as much as 20 percent and 60 percent, respectively, with temperature differ- 
ences as large as those between the frozen shock temperature and the temperature 
at the [0] [CO] maximum ( =s 200 K) . For this reason, the temperature predicted by 
the kinetic model, T m in table I, was used to calculate [CO] and [[002]. The 
justification for the use of T m is discussed in detail in reference 5. 
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COMPARISON BETWEEN EXPERIMENTAL RESULTS AND KINETIC MODELS 


Propane decomposition and oxidation behind shock waves have been investi- 
gated in previous shock-tube experiments (refs. 6 to 8). One of these investiga- 
tions (ref. 6) concentrated on the ignition characteristics of propane. Another 
(ref. 8) concentrated on the decomposition of propane behind a reflected shock 
wave and suggested a sequence of reaction steps and rate constants which analyti- 
cally predicted the experimental measurements of the investigation. Using this 
information, a reaction scheme was assembled which could be studied parametri- 
cally to determine the influence of individual steps in achieving a good compari- 
son with the experimental measurements of this paper as well as with ignition 
delay time data from reference 6. The preliminary kinetic mechanism with the 
associated rate coefficients is listed in table II. The rate coefficients were 
obtained from the literature whenever possible. (See refs. 4, 5, 8, 1 3 » 14, and 
15.) The rate coefficient for reaction (1) was estimated by using the value 
reported by Lifshitz and Frenklach (ref. 8) at 1200 K and assuming an activa- 
tion energy equal to the bond dissociation energy. The rate coefficients for 
reactions (6) and (7) were estimated using Semenov's rule (ref. 16) to deter- 
mine the activation energy and the Arrhenius factor was arbitrarily set at 
5.00 x 1 o 1 3 cm3/mol-sec. For reaction ( 1 4 ) , a zero activation energy was 
assumed since the reaction is highly exothermic, and the rate coefficient was 
set at 10^3 cm3/mol-sec. 

The analytical model used in the parametric study involved one-dimensional 
flow behind an incident shock wave with finite-rate chemical reactions and 
boundary-layer growth. The basic computer program is described in reference 12; 
however, a modified form of this computer program described in reference 17 was 
used for most of the calculations because of its faster computational times for 
systems containing many species and reactions. The thermochemical data of 
JANAF (ref. 18) were used for the species involved in the mechanism; however, 
when JANAF data were not available, the data of reference 19 were used. 

The parametric study was performed to obtain a good comparison between the 
present experimental results for the [o] [CO] product and the product calculated 
by the analytical model. Only the rate constants involving propane (reac- 
tions (1) to (14)) were manipulated within the suspected error limits. At the 
same time the additional restraint was imposed that the same rate mechanism must 
reproduce reasonably well the ignition delay time data presented in reference 6. 
The ignition delay time calculations were made with the same computer programs 
(refs. 12 and 17) operating in a manner which would simulate conditions behind a 
reflected shock wave, the condition under which the ignition delay time data of 
reference 6 were taken. 

The first objective of the parametric study was to identify reactions which 
significantly influence the predicted [6] [CO] profiles, reaction times t m , and 
ignition delay times t^. To do this, the rate coefficients for reactions (1) 
and (14) were varied by a factor of 100, for reactions (2), (3), (10), (11), 

(12), and (13) by a factor of 10, and for reactions (4) to (9) by a factor of 4. 
The predicted [0] [CO] profiles, reaction times, and ignition delay times were 
found to be most sensitive to the rates of the following reactions: 
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( 1 ) 

c 3 h 8 + 

ch 3 + 

c 2 h 5 

( 2 ) 

ch 3 + 

c 3 h 8 - 

CH 4 + n-C 3 Hy 

(3) 

ch 3 + 

c 3 h 8 + 

CH|j + i-C 3 H 7 

( 10 ) 

n-C 3 Hy 


+ ch 3 

( 12 ) 

n-C 3 Hy 

- c 3 h 6 

+ H 

(13) 

i-C 3 Hy 

- c 3 H 6 

+ H 


The predicted [ 6 ] [CCQ profiles and reaction times, however, were influenced most 
by the rate of reaction ( 1 ), whereas the predicted ignition delay times were 
influenced most by the rates of reactions (1), (2), and (3). The reaction mech- 
anism and rate coefficients listed in table II predicted [_0] [CO] values that 
were about 50 percent larger than the measured values and reaction times that 
were about a factor of 2 smaller. The predicted ignition delay times were from 
10 to 50 percent smaller than the measured values. 

In the final phase of the analytical study, an attempt was made to find a 
set of rate coefficients for the six reactions that would reproduce as well as 
possible the experimental data without resorting to unreasonable rate coeffi- 
cients. Initially, the rate coefficient for reaction (1) was adjusted to pro- 
vide reasonable agreement between the predicted and measured [b] (jCO] values 
while the rate coefficients for reactions ( 2 ), ( 3 ), ( 10 ), ( 12 ), and ( 13 ) were 
fixed. It was found that best agreement was obtained when the decomposition of 
propane was written as a bimolecular process 

( 1 ) M + C 3 H 8 + CH 3 + C 2 H 5 + M 

with the rate coefficient, k-j □ 5.00 * 10^5 exp(-32 713/T) cm3/mol-sec. The 
need to describe the decomposition of propane as a second-order process with an 
activation energy less than the bond dissociation energy suggests that the decom- 
position reaction is in the fall-off region of a unimolecular decomposition pro- 
cess. Similar behavior has been observed for methane (ref. 20) and methanol 
(ref. 21). The measured and predicted reaction times were also in good 
agreement . 

After the rate coefficient for reaction (1) had been established, an 
attempt was made to find a set of rate coefficients for the other reactions that 
reproduced as well as possible the measured ignition times. With the adjusted 
rate coefficient for reaction ( 1 ) , the predicted ignition times were substan- 
tially larger (a factor of 2 or greater) than the measured values. The ignition 
times could be decreased in four ways: by increasing the rates of reactions ( 2 ) 

and ( 3 ), by decreasing the rate of reaction ( 10 ), by increasing the rates of 
reactions (12) and (13), or by a combination of the above. Each of the four 
ways was attempted; however, a set of rate coefficients could not be found that 
were considered reasonable and that did not significantly influence the pre- 
dicted [0] [CO] values and reaction times. Therefore, it appeared reasonable to 
reexamine the proposed reaction scheme to determine the cause of the long pre- 
dicted ignition delay times. After examining several computer runs, it became 
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obvious that the long ignition delays were caused by reaction (15) which tied up 
the hydrogen radical, and since no removal reaction was present other than decom- 
position, there was a large buildup of the HO 2 radical. To remedy this, two 
reactions were added to the preliminary scheme (see table III): 

(58) H + H0 2 OH + OH 

(59) CO + H0 2 -► OH + C0 2 

The rate coefficient for reaction (58) was set at (ref. 13) * 

2.5 x 10 ^ exp(-950/T) cm3/mol-sec v 

An experimental rate coefficient for reaction (59) could not be found; however, 
a recent study (ref. 22) of reaction (25), CO + OH -*■ C0 2 + H, suggested that 
the rate coefficient for reaction (59) must be significantly smaller than the 
rate coefficient for reaction (25). Therefore, the rate coefficient for reac- 
tion (59) was set at 10^ cm^/mol-sec which is about a factor of 10 smaller 
than the rate coefficient for reaction (25) over the temperature range of this 
study . 

By incorporating these reactions into the reaction scheme, the predicted 
ignition delay times were substantially smaller, especially at lower tempera- 
tures but still about 50 percent higher than the measured values. The predicted 
[ 0 ] [CO] values and reaction times were not affected by these reactions. Better 
agreement between the measured and predicted- ignition times was obtained by 
adjusting the rate coefficient of reactions (2) and (3) to 

2.00 * 10 ^ exp(-5l84/T) cm3/mol-sec 

No attempt was made to vary the rate coefficients between reactions (2) and (3). 

This rate coefficient is about a factor of 5 larger than the value reported by 
Lifshitz and Frenklach (ref. 8) but is not an unreasonable rate coefficient for 
this type of reaction. A comparison of the predicted and measured [ 0 ] [CO] val- 
ues is given in figure 7, and of the ignition delay times in table IV. Figure 8 
shows the comparison between the predicted and measured concentrations of CO and 
C0 2 over the temperature range of this study. Figure 9 shows the comparison of 
predicted and measured reaction times t m . In all the comparisons, the values 
predicted using the final kinetic mechanism agreed well with the experimental 
data. 

/ 

The final kinetic mechanism for propane is listed in table III. The entire 
mechanism consists of 59 kinetic steps. To reduce the number of kinetic steps v 

for the propane mechanism, it was found that the following reactions can be 
omitted from the mechanism without changing the computed values: 

(36) C 2 H 4 + C 2 H 4 -> c 2 h 5 + c 2 h 3 

(35) C 2 H 4 + M C 2 H 3 + H 

(46) C 2 H 2 + 0 2 -*• HCO + HCO 
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(45) C 2 H 2 + M C 2 H + H 


These reactions were carried over from the kinetic mechanisms for acetylene and 
ethylene from reference 5. 


CONCLUDING REMARKS 

The stoichiometric combustion of propane behind incident shock waves was 
studied experimentally and analytically over a temperature range from 1700 K to 
2600 K and a pressure range from 1.2 to 1.9 atm. Measurements of the concentra- 
tions of carbon monoxide [CO] and carbon dioxide [C0 2 ] and the product of the 
oxygen atom and carbon monoxide concentrations were made after passage of the 
incident shock wave. Analytical comparisons with the experimental data were 
made using a computer model of a flowing, reacting gas behind incident shock 
waves. In order to use the analytical model, a kinetic mechanism was assembled 
based on the assumption that propane decomposes rapidly to fragments with lower 
molecular weights followed by the oxidation of these fragments. A parametric 
study identifying and adjusting the rate coefficients of the reaction steps most 
sensitive in predicting the [o] [CO] concentration product, reaction times, and 
ignition delay times led to a kinetic mechanism which predicted well all present 
experimental results and experimental ignition delay times of a previous shock- 
tube study. The agreement between the computed and measured ignition delay 
times was also good. The kinetic mechanism, and associated rate coefficients, 
consists of 59 steps. Almost all the rate coefficients were taken from the 
literature. 

The assumption that the parent hydrocarbon, propane in this case, decom- 
poses rapidly to fragments with lower molecular weights followed by the oxida- 
tion of these fragments during later stages of the combustion process appears 
to be a valid approach. Since recent literature suggests that high-molecular- 
weight hydrocarbons such as iso-octane, toluene, and cyclohexane tend to follow 
this type of reaction sequence, there is a need to determine whether similar 
kinetic models can be assembled for these fuels. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
May 18, 1977 
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TABLE I.- EXPERIMENTAL RESULTS 


(a) At observation station (T) 


T, 

P, 

[0] [CO] , 

K 

atm 

( mol / cm 3)2 

2525 

1 .81 

9.70 X 10-15 

2522 

1.81 

9.70 

2510 

1.80 

10.2 

2510 

1 .80 

11.2 

2505 

1.80 

96.7 

2470 

1.77 

10.4 

2422 

1 .74 

11.2 

2296 

1 .62 

— 

2277 

1 .60 

— 

2247 

1 .59 

8.65 

2247 

1 .59 

8.59 

2235 

1 .56 

8.76 

2073 

1 .44 

6.82 

1987 

1 .37 

7.44 

1887 

1 .28 

5.64 

1887 

1.28 

5.60 

1880 

1 .28 

6.53 

1874 

1 .27 

5.92 

1870 

1 .27 

5.57 

1862 

1.26 

6.38 


^m» 

[c° 2 ]. 

[CO], 

ysec 

mol/cnP 

mol/cmP 

4.0 

_ _ _ _ 

1.69 x 10- 

4.0 

— 

1 .70 

4.0 

— 

1.71 

4.0 

4.40 X icr 8 

— 

4.0 

4.59 

— 

4.0 

4.30 


4.0 

4.62 

— 

— 

4.58 

— 

— 

— 

1.79 

5.0 

4.89 

— 

5.0 

4.46 


6.0 

— 

1 .76 

9.0 

— 

1 .67 

9.5 

— 

1 .83 

10.0 

— 

1 .34 

12.5 

3.76 



12.0 

— 

1 .52 

11.0 

3.85 

— 

11 .0 

3.83 

— 

10.0 

— 

1 .55 



2710 

2709 

2697 

2697 

2692 

2660 

2616 

2503 

2486 

2459 

2459 

2448 

2309 

2237 

2157 

2157 

2152 

2146 

2143 

2137 
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TABLE I.- Concluded 


Cb) At observation station (D 


T, 

K 

P, 

atm 

[o] [col , 

(mol/cm3) 2 

ysec 

2615 

1.90 

10.5 x 10-15 

3.50 

2570 

1.86 

9.18 

2.80 

2480 

1.77 

10.1 

4.00 

2415 

1.71 

8.99 

3.50 

2300 

2300 

1.61 

1.61 

8.72 

8.21 

5.50 


2255 

1.57 

8.31 

5.50 

2170 

1.50 

8.50 

7.50 

2180 

1.51 

7.95 

7.00 

2005 

1.37 

8.38 

10.0 

1190 

1895 

1805 

1.36 

1.28 

1.21 

7.86 

6.82 

6.31 

9.00 

12.0 



[“>£ • , 
mol/cnP 

[CO] 

mol/cnw 

I'm > 
K 

4.45 x 10“ 8 

1.48 x 10"7 

2793 

4.20 

1.49 

2751 

4.76 

1 .74 

2670 

4.24 

1 .55 

26 11 

4.15 

1.72 

2507 

3-78 

1.61 

2507 

4.26 

1.72 

2467 

4.07 

1.77 

2392 

3-68 

1.72 

2401 

4.07 

1.48 

2254 

4.12 

1.65 

2240 

3.65 

1 .63 

2165 

3.88 

1.74 

2095 



TABLE II.- PRELIMINARY KINETIC MECHANISM FOR PROPANE COMBUSTION 


Reaction 

Rate coefficient, k 
(a) 

Ref. 

( 1 ) c 3 h 8 -*■ ch 3 + c 2 h 5 

5.00 x i 0 15 exp(-43 785/T) 

Estimated 

(2) CH 3 + C 3 Hg -*■ CH 4 + n-C 3 H 7 

3-55 x i 0 12 exp(-5l84/T) 

8 

(3) CH 3 + C 3 H 8 CH 4 + i-C 3 Hy 

3-55 x i 0 12 exp(-5l84/T) 

8 

(4) H + C 3 H 8 •> H 2 + n-C 3 H 7 

6.30 x 10 1 3 eX p(-4026/T) 

8 

(5) H + C 3 H 8 -*■ H 2 + i-C 3 H 7 

6.30 x 10 1 3 eX p(-4026/T) 

8 

(6) 0 + C 3 H 8 -> OH + n-C 3 H 7 

5.00 x 10 1 3 ex p(-5033/T) 

Estimated 

(7) 0 + C 3 H 8 ■> OH + i-C 3 H 7 

5-00 x 1 0 1 3 ex p(-5033/T) 

Estimated 

(8) OH + C 3 H 8 -* H 2 0 + n-C 3 H 7 

1.60 x 10 14 exp(-1580/T) 

15 

(9) OH + C 3 H 8 -> H 2 0 + i-C 3 H 7 

1.60 x 10 14 exp(-1580/T) 

15 

(10) n-C 3 H 7 ■> C 2 H i, + CH 3 

4.00 x 10 1 3 exp (-16 658/T) 

14 

(11) i-C 3 H 7 C 2 H|, + CH 3 

1.00 x 1012 exp (-17 363/T) 

8 

(12) n-C 3 H 7 -> C 3 H 6 + H 

6.30 x 10 1 3 exp (-19 124/T) 

14 

(13) i-C 3 H 7 + C 3 H 6 + H 

2.00 x 10 121 exp (-20 785/T) 

14 

(14) 0 + C 3 H 6 -* CH 2 0 + C 2 H 4 

1.00 x 1013 

Estimated 

(15) M + H + 0 2 -> H0 2 + M 

1.50 x 10 1 ^ eX p(500/T) 

13 

(16) M + CH 4 ->■ CH 3 + H + M 

4.00 x 1017 exp (-44 500/T) 

4 

(17) CH 4 + H -»■ CH 3 + H 2 

1.26 x 10 14 exp(-5989/T) 

4 

(18) CH 4 + OH -> CH 3 + H 2 0 

3.00 x 10 1 3 ex p(-3020/T) 

4 

(19) CH 4 + 0 CH 3 + OH 

2.00 x 1 0 1 3 eX p(-4640/T) 

4 


a The units for k are sec - 1 for unimoleeular reactions, cm3/mol-sec 
for bimolecular reactions, and cm^/mol^-sec for termolecular reactions. 
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TABLE II.- Continued 


Reaction 

Rate coefficient, k 

Ref. 


(a) 


(20) CH 3 + 0 2 -»■ CH 2 0 + OH 

1.70 x 10 12 exp(-7045/T) 

4 

(21 ) . CH 3 + 0 CH 2 0 + H 

1.30 x 10 14 exp(-1006/T) 

4 

(22) CH 2 0 + OH + HCO + H 2 0 

2.30 x ioU 

4 > 5 

(23) HCO + OH -*■ CO + H 2 0 

1.00 x 10 11 * 

5 

(24) M + HCO -*■ H + CO + M 

5.00 x 10 14 exp(-9562/T) 

5 

(25) CO + OH -»• C0 2 + H 

4.00 x 10 12 exp(-4026/T) 

4 , 5 

(26) CO + 0 + M -*■ C0 2 + M 

6.00 x loU 

5 

(27) 0 + H 2 ->■ OH + H 

2.07 x 10 14 exp(-6920/T) 

4, 5 

(28) H + 0 2 -*• OH + 0 

1.22 x io 1 7t-°- 91 exp(-8369/T) 

4 , 5 

(29) H 2 + OH -»■ H 2 0 + H 

5.20 x loU exp(-3271/T) 

4, 5 

(30) OH + OH -► H 2 0 + 0 

5.50 x ioU exp(-3523/T) 

4 > 5 

( 31 ) h + h + m-*h 2 + m 

1.00 x lO 1 ^ 

5 

(32) H + OH + M -* H 2 0 + M 

8.40 x io 21 T - 2 - 0 

4, 5 

(33) 0 2 + H 2 ■> OH + OH 

1.70 x 10 U exp(-24 230/T) 

4 

(34) M + 0 2 ->-0 + 0 + M 

2.55 x ioISt - 1 * 0 exp (-59 386 /T) 

4 , 5 

(35) M + C 2 H 4 -*■ C 2 H 3 + H + M 

1.00 x 10 14 exp(-54 857/T) 

5 

(36) C 2 H 4 + C 2 H 4 -»• C 2 H 5 + C 2 H 3 

5.00 x 10 14 exp(-32 562/T) 

5 

( 37 ) C 2 H 5 -*■ C 2 H 4 + H 

3.16 x loU exp (-20 483/T) 

5 

(38) M + C 2 H 3 -*■ C 2 H 2 + H + M 

3.00 x 10 16 exp (-20 382/T) 

5 


a The units for k are sec“^ for unimolecular reactions, cm3/mol-sec 
for bimolecular reactions, and cm^/mol 2 -sec for termolecular reactions. 


15 


TABLE II.- Concluded 


Reaction 

Rate coefficient, k 

Ref 

— — _ 

(a) 


(39) H + C 2 H 4 -*■ C 2 H 3 + H 2 

1.10 X 10 14 exp(-4278/T) 

5 

(40) OH + C 2 H 4 -> C 2 H 3 + H 2 0 

1.00 x 10 14 exp(-176l/T) 

5 

(41 ) 0 + C 2 H 4 -»• CH 2 0 + CH 2 

2.50 x lO^ exp(-25l6/T) 

5 

(42) 0 + C 2 H 4 -*• CH 3 + HCO 

2.26 x 10^ exp(-1359/T) 

5 

(43) CH 2 0 + H -*• HCO + H 2 

2.00 x 10^ 

5 

(44) CH 2 0 + 0 -> HCO + OH 

2.00 x 10^ 

5 

(45) M + C 2 H 2 -»■ C 2 H + H + M 

1.00 x 10 14 exp (-57 373/T) 

5 

(46) C 2 H 2 + 0 2 •> HCO + HCO 

4.00 x 10 12 exp(-l4 092/T) 

5 

(47) H + C 2 H 2 -»■ C 2 H + H 2 

2.00 x TO 11 * exp<-9562/T) 

5 

(48) OH + C 2 H 2 -*■ C 2 H + H 2 0 

6.00 x 10 12 exp(-3523/T) 

5 

(49) 0 + C 2 H 2 -»■ C 2 H + OH 

3-20 x 10 15 T-°-6 exp(-8556/T) 

5 

(50) 0 + C 2 H 2 ■+ CH 2 + CO 

5.20 x 10 ^ exp(-l862/T) 

5 

(51) C 2 H + 0 2 HCO + CO 

1.00 x 10 ^ exp(-3523/T) 

5 

(52) C 2 H + 0 -*• CO + CH 

5.00 x 10^3 

5 

(53) CH 2 + 0 2 -► HCO + OH 

1.00 x 10 14 exp(-l862/T) 

5 

(54) CH + 0 2 -* CO + OH 

1.35 x ioDt 0 * 6 ^ exp(-12 934/T) 

5 

(55) CH + 0 2 -► HCO + 0 

1.00 x 10*13 

5 

(56) HCO + H ■*> CO + H 2 

1.00 x 10 14 

5 

(57) HCO + 0 -*■ CO + OH 

1.26 x I0 11j 

5 


a The units for k are sec - ^ for unimolecular reactions, 
cm3/mol-sec for bimolecular reactions, and cm^/mol 2 -sec for termo- 
lecular reactions. 



TABLE III.- FINAL KINETIC MECHANISM FOR PROPANE COMBUSTION 


Reaction 

Rate coefficient, k 
(a) 

Ref. 

<> 

( 1 ) M + C 3 H 8 -*• C 2 H 5 + CH 3 + M 

5.00 x 10*15 exp(-32 713/T) 

Present study 

(2) CH 3 + C 3 H 8 CH 4 + n-C 3 H 7 

2.00 x lO*^ exp(-5l84)/T) 

Present study 

(3) CH 3 + C 3 H 8 ■> CH 4 + i-C 3 H 7 

2.00 x lO^ exp(-5l84)/T) 

Present study 

(4) H + C 3 H 8 H 2 + n-C 3 H 7 

6.30 x lO*^ exp(-4026/T) 

8 

(5) H + C 3 H 8 ■> H 2 + i-C 3 H 7 

6.30 x i 0 13 exp(-4026/T) 

8 

( 6 ) 0 + C 3 H 8 -*■ OH + n-C 8 H 7 

5.00 x 10 U exp(-5033/T) 

Estimated 

(7) 0 + C 3 H 8 -*■ OH + i-C 3 H 7 

5.00 x lO^ exp(-5033/T) 

Estimated 

( 8 ) OH + C 3 H 8 •> H 2 0 + n-C 3 H 7 

1.60 x lO 1 * 1 exp(-1580/T) 

15 

(9) OH + C 3 H 8 -*• H 2 0 + i-C 3 H 7 

1 .60 x lO 14 exp(-1580/T) 

15 

(10) n-C 3 H 7 -*■ C 2 H 4 + CH 3 

4.00 x lO^ exp (-16 658 /T) 

14 

(11) i-C 3 H 7 C 2 H 4 + CH 3 

1.00 x 10 12 exp (-17 363/T) 

8 

(12) n-C 3 H 7 ->• C 3 Hg + H 

6.30 x 10*>3 exp (-19 124/T) 

14 

(13) i-C 3 H 7 C 3 Hg + H 

2.00 x 10 14 exp (-20 785/T) 

14 

(14) 0 + C 3 H 6 CH 2 0 + C 2 H 4 

1.00 x 10 U 

Estimated 

(15) M + H + 0 2 -* H0 2 + M 

1.50 x 10*15 eX p(500/T) 

13 

(16) M + CH 4 -*■ CH 3 + H + M 

4.00 x 1017 exp(-44 500/T) 

4 

(17) CH 4 + H -»■ CH 3 + H 2 

1.26 x 10 14 exp(-5989/T) 

4 

(18) CH 4 + OH -*■ CH 3 + H 2 0 

■3-00 x 10 ^ exp(-3020/T) 

4 

(19) CH 4 + 0 ♦ CH 3 + OH 

2.00 x 10 ^ exp(-4640/T) 

4 


a The units for k are sec"1 for unimolecular reactions, cm3/mol-sec 
for bimolecular reactions, and cm^/mol^-sec for termolecular reactions. 


TABLE III.- Continued 


Reaction 

Rate coefficient, k 
(a) 

Ref. 

(20) CH 3 + 0 2 -»■ CH 2 0 + OH 

1.70 x 10 12 exp(-7045/T) 

4 

(21 ) CH 3 + 0 ■* CH 2 0 + H 

1-30 x 10 11 * exp(-1006/T) 

4 

(22) CH 2 0 + OH -*• HCO + H 2 0 

2.30 x lO^ 

4, 5 

(23) HCO + OH ->■ CO + H 2 0 

^r 

0 

X 

0 

0 

5 

(24) M + HCO H + CO + M 

5.00 X lO 1 * 1 exp(-9562/T) 

5 

(25) CO + OH + C0 2 + H 

4.00 x 10 12 exp(-4026/T) 

it, 5 

(26) CO + 0 + M -»■ C0 2 + M 

6.00 x 10 1 3 

5 

(27) 0 + H 2 ->■ OH + H 

2.07 x lO 14 exp(-6920/T) 

4, 5 

(28) H + 0 2 -*■ OH + 0 

1.22 x io17t- 0.91 exp(-8369/T) 

5 

(29) H 2 + OH -*• H 2 0 + H 

5.20 x 10 "I 3 exp(-3271/T) 

4, 5 

(30) OH + OH ^ H 2 0 + 0 

5.50 x lO^ exp(-3523/T) 

5 

(31) h + h + m-»-h 2 + m 

LD 

O 

X 

O 

O 

5 

(32) H + OH + M -*■ H 2 0 + M 

8.40 x io 21 T ~ 2 - 0 

4 , 5 

(33) 0 2 + H 2 -*■ OH + OH 

1.70 x 10^3 exp (-24 230/T) 

4 

(34) M+0 2 -*-0 + 0+ M 

2.55 x io i 8 T-1 *° exp (-59 386/T) 

J t, 5 

(35) M + C 2 H 4 + C 2 H 3 + H + M 

1.00 x lO 14 exp(-54 857/T) 

5 

(36) C 2 Hj| + C 2 H 4 c 2 h 5 + c 2 h 3 

5.00 x lO 12 * exp(-32 562/T) 

5 

(37) C 2 H 5 -► C 2 H 4 + H 

3. 16 x 1013 exp (-20 483/T) 

5 

(38) M + C 2 H 3 ■* C 2 H 2 + H + M 

3-00 x 10 16 exp (-20 382/T) 

5 


a The units for k are sec -1 for unimolecular reactions, cm3/mol-sec 
for bimolecular reactions, and cm°/mol 2 -sec for termolecular reactions. 
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TABLE III.- Concluded 


Reaction 

Rate coefficient, k 
(a) 

Ref. 

(39) H + C 2 H 4 -»• C 2 H 3 + H 2 

1.10 x 10 14 exp(-4278/T) 

5 

(40) OH + C 2 H 4 C 2 H 3 + H 2 0 

1 .00 x lO 14 exp(-176l/T) 

5 

(41) 0 + C 2 Hh -+ CH 2 0 + CH 2 

2.50 x lO^ exp(-25l6/T) 

5 

(42) 0 + C 2 H 4 CH 3 + HCO 

2.26 x lO^ exp(-1 359/T) 

5 

(43) CH 2 0 + H ->• HCO + H 2 

2.00 x lO^ 

5 

(44) CH 2 0 + 0 -»• HCO + OH 

2.00 x lO^ 

5 

(45) M + C 2 H 2 C 2 H + H + M 

1.00 x lO 14 exp(-57 373/T) 

5 

(46) C 2 H 2 + 0 2 HCO + HCO 

4.00 x 10 12 exp (-14 092/T) 

5 

(47) H + C 2 H 2 -► C 2 H + H 2 

2.00 x 10 14 exp(-9562/T) 

5 

(48) OH + C 2 H 2 -v C 2 H + H 2 0 

6.00 x 10 12 exp(-3523/T) 

5 

(49) 0 + C 2 H 2 ->■ C 2 H + OH 

3.20 x io15t- 0.6 exp(-8556/T) 

5 

(50 ) 0 + C 2 H 2 -> CH 2 + CO 

5.20 x lO^ exp(-l862/T) 

5 

(51) C 2 H + 0 2 ->• HCO + CO 

1.00 x 10^ exp(-3523/T) 

5 

(52) C 2 H + 0 + CO + CH 

5.00 x lO^ 

5 

(53) CH 2 + 0 2 + HCO + OH 

1.00 x 1 0 1 ^ exp(-l862/T) 

5 

(54) CH + 0 2 ■+• CO + OH 

1.35 x ioIIt 0 - 6 ^ exp (-12 934/T) 

5 

(55) CH + 0 2 ->• HCO + 0 

1 .00 x lO^ 

5 

(56) HCO + H ->■ CO + H 2 

1.00 x lO 14 

5 

(57) HCO + 0 -► CO + OH 

1.26 x lO 14 

5 

(58) H + H0 2 -> OH + OH 

2.50 x 10 14 exp(-950/T) 

19 

(59) CO + H0 2 + OH + C0 2 

1.00 x 10 10 

Estimated 


a The units for k are sec - ^ for unimolecular reactions, cm3/mol-sec 
for bimolecular reactions, and cm 6 /mol 2 -sec for termolecular reactions. 
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TABLE IV.- COMPARISON BETWEEN MEASURED AND CALCULATED 


IGNITION TIMES 


Gas composition, 
percent mole 

T, K 

p, atm 

V 

ysec 

c 3 h 8 

°2 

Ar 

Measured 

Calculated 

1.60 

8.00 

90.40 

1327 

8.59 

282 

320 

1 .60 

8.00 

90.40 

1352 

8.45 

185 

220 

.80 

8.00 

VO 

[\> 

o 

1423 

7.66 

70 

65 

1.60 

4.00 

94.40 

1447 

8.10 

225 

300 

.41 

4.10 

95.49 

1492 

7.67 

77 

60 

1 .60 

8.00 

90.40 

1495 

2.62 

68 

75 

.84 

2.10 

97.06 

1539 

8.32 

230 

230 

3-85 

19-23 

76.92 

1587 

14.07 

12 

10 


a From reference 6. 
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Figure 1 .- Photograph of shock tube indicating location of 
test observation stations. 
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Schematic representation of instrumentation arrangement 
at observation stations. 
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(b) Observation station (?) . 
Figure 2.- Concluded. 
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(a) Variation of CO emission voltage at 5.0 urn with temperature. 



(b) Correlation of molar concentration with temperature and emission voltage. 


Figure 4.- Calibration of detection system for obtaining CO molar 

concentration , [CO] . 
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Figure 6.- Typical emission profiles at 4.3 pm, 5.0 pm, and 3700 A. 
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Figure 7 


- Comparison of predicted and measured [o] [CO] values from 
stoichiometric propane combustion. 
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Figure 8.- Comparison of predicted and measured molar concentrations 
of CO and CO 2 during stoichiometric combustion of propane. 
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